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Abstract We investigate the thermoelectric properties of PbTe doped with a small 
concentration x of Tl impurities acting as acceptors and described by Anderson im- 
purities with negative on-site (effective) interaction. The resulting charge Kondo 
effect naturally accounts for a number of the low temperature anomalies in this 
&> | system, including the unusual doping dependence of the carrier concentration, the 



Fermi level pinning and the self-compensation effect. The Kondo anomalies in the 

Slow temperature resistivity at temperatures T < 10 K and the x-dependence of the 
residual resistivity are also in good agreement with experiment. Our model also cap- 
""^J | tures the qualitative aspects of the thermopower at higher temperatures T > 300 K 

Ch ■ for high dopings (x > 0.6%) where transport is expected to be largely dominated by 

carriers in the heavy hole band of PbTe. 



>• ■ 1 Introduction 

f^i | The ever growing demand for energy has increased interest in materials with high 

C*~) ■ thermoelectric (TE) efficiency. A remarkable feature of TE devices is that they di- 

tH- J rectly convert heat into electrical energy or use electricity to pump heat by circu- 

C^ ■ lating the electron fluid between a hot and a cold reservoir (jT] |2). The absence of 

C*~) | mechanical parts makes TE devices very reliable but their low efficiency restricts 

applications to rather specialized fields. Today, they are mainly used in situations 
where reliability is the most important factor, e.g., for electricity generation in re- 
Kjf . mote regions or for vibrationless cooling of sensitive devices. 
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The efficiency of TE materials is related to their figure-of-merit, ZT = P T /k, 
where P is the power factor and K is the thermal conductivity. Current applications 
are based on semiconducting materials that have a state-of-the-art ZT of around 
one 031, which is not sufficient for wide-spread use. The production of cheap TE 
materials with a room temperature value of ZT > 2 would be a major technological 
breakthrough. 

The efforts to increase ZT follow two general strategies. The one tries to reduce K 
and increase the efficiency by minimizing the thermal losses through the device. The 
other tries to increase P = S 2 (7, were S(T) is the Seebeck coefficient and <j(T) the 
electrical conductivity. Recent advances that achieved ZT > 2 are mostly due to the 
reduction of K by nanostructuring, i.e., by using multilayered materials (for review 
see (4))andnanocomposites||5]. The enhancement of ZT via the power factor is less 
straightforward, because an increase of S(T) is usually accompanied by a reduction 
of a and vice versa |6|. The reason is that o(T) is mainly determined by the value 
of the density of states ^K(o) and the transport relaxation time t(g>) at the chemical 
potential ji, while S(T) depends sensitively on the symmetry of jV(a>) and t(<b) 
around fi. Thus, one way to a larger power factor is by increasing the asymmetry of 
JfifQ) in the proximity of \i by, say, confining the electron gas to low dimensions 
Q or by introducing resonant states within the Fermi window by doping [8 1. This 
strategy could work well for materials in which the charge carriers behave as non- 
interacting Fermions and the energy dependence of x{(0) is negligibly small J9j- 

The other route to large P is via electron correlation which makes t(co) strongly 
energy dependent, as in heavy fermions J6] [10], e.g., YbAl3, transition metal ox- 
ides, e.g. Na x Co02 iTTTlfETl . or in correlated semiconductors, such as FeSi lfi"3"1[l4l . 
However, these materials have large K such that ZT is too small for applications, 
except, perhaps, at very low temperatures. At high temperatures, the correlation ef- 
fects are usually less pronounced and the appropriate materials are not easy to find. 
Considering the complexity of the problem, it is clear that the search for better ther- 
moelectrics could benefit considerably from an understanding of basic microscopic 
mechanisms that govern the thermoelectric response of a given material. 

One of the most interesting thermoelectric materials is Pbi_ A TLTe, where x is 
the concentration of Tl impurities 1131 [T6l [T71 . The parent compound PbTe is a nar- 
row gap semiconductor with a gap of 190meV and is itself a good thermoelectric 
material with a thermopower of around 300/iV/K at 800 K. Upon doping with Tl 
impurities, which act as acceptors lfT5] [T6) . a large thermopower persists at high 
temperature, while a number of anomalies appear in the low temperature proper- 
ties. Remarkable, these depend sensitively on the Tl concentration with a quali- 
tatively different behavior below and above a critical concentration x c ~ 0.3 at% 
Tl. For example, Pbi_ x Tl v Te becomes superconducting for x > x c with a transition 
temperature T c (x) increasing linearly with x and reaching 1.5K at x = 1.5 at% Tl 
|[T8l[T9ll . This is surprisingly high given the low carrier density of 10 20 holes/cm 3 . 
Measurements of the Hall number p# = 1 /Rne |T9l indicate that the number of 
holes grows linearly with x for x < x c , whereas for x > x c the number of holes re- 
mains almost constant: the system exhibits "self-compensation" and the chemical 
potential is pinned to a value ji ~ jj,* = 220 meV lfl9l l20l I2TI . Transport measure- 
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merits show anomalous behavior at low temperatures: while forx < x c , the residual 
resistance po is very small and almost constant J22), forx > x c , po increases approx- 
imately linearly with x. Forx < x c , the resistivity, p(T) exhibits a positive slope at 
low temperature El . while for x > x c , the slope is negative and a Kondo like im- 
purity contribution pj mp (r) is observed for T < 10K. The origin of this anomaly is 
not due to magnetic impurities, since the susceptibility is diamagnetic |23l . Finally, 
at low Tl concentrations, S(T) shows a sign-change before growing to large positive 
values at room temperatures El . 

Several models have been proposed to explain the anomalous properties of 
Pbi-^TljcTe. The simplest assumes that doping with Tl gives rise to non-interacting 
resonant levels close to the top of the p-valence band of PbTe |fl31 . Density func- 
tional theory (DFT) calculations confirm the presence of such states with Tl s- 
character El . In the static mixed valence model El [TBI , Tl impurities, known 
to be valence skippers in compounds, are assumed to dissociate into energetically 
close Tl 1+ (6s 2 p 1 ) and Tl 3+ (6s°p 3 ) ions, while the Tl 2+ (6s l p 2 ) configuration lies 
higher in energy (consistent with DFT calculations [26]). With a strong electron- 
phonon interaction such a model results in negative on-site U. This model provides 
a natural explanation for the observed self-compensation and the diamagnetic be- 
havior of Pbi_^Tl^Te El . In addition, it provides a mechanism for the onset of 
superconductivity. Dynamic fluctuations between the Tl 1+ and Tl 3+ valence states 
results in the negative U Anderson model El which supports a charge Kondo effect 
Il28ll29l and has been proposed to explain the anomalous properties of Pbi_ r Tl A Te 
EJ [HI [22] [30] EQ E2 ED, including the superconductivity for x > x c EQ ED . 
A coupling of the Tl ions to the lattice has also been considered l34l [35] . Recent 
ARPES data on a 0.5% sample J36l are also consistent with a negative U Anderson 
model. 

As regards the thermoelectric properties of interest for applications, increasing 
the Tl concentration beyond 1 at% gives rise to a room-temperature ZT which is 
surprisingly large for a bulk material |fT7| . The data show that the thermopower 
S(T) increases almost linearly at low temperatures El and saturates above 400 K 
at rather high values 0171 . The fact that the enhancement of ZT occurs at high tem- 
peratures and that it is mainly due to an increased power factor makes this material 
particularly interesting not just from the practical but also from the theoretical point 
of view. 

In this paper, we use the numerical renormalization group approach (NRG) l37l 
to calculate the anomalous frequency and temperature dependent transport time, 
t(co,T), of electrons scattering from Tl impurities described as negative-C/ Ander- 
son impurities. We show that this model explains the unusual concentration and 
temperature dependent properties of Pbi_ v Tl x Te. We take full account of charge 
neutrality and include the realistic band structure of PbTe |9|. We find that both cor- 
relation and band effects are important in accounting for the experimental data. We 
shall also show, that our single-band Anderson impurity model qualitatively repro- 
duces the measured behaviour of the thermopower over a wide temperature range 
in the regime where transport is mainly due to carriers in a single heavy hole like 
sub-band of the PbTe valence band. At lower dopings, the population of two valence 
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sub-bands of PbTe requires a further generalization of the model in order to describe 
the thermopower also in this regime (see discussion in Sec. |5J, This, however, is be- 
yond the scope of the present paper. 



2 Negative-^/ Anderson Model for Tl impurities in PbTe 

We consider n Tl impurities in a PbTe crystal with N Pb sites described by the 
Hamiltonian H = i/band + #imp + #hyb, where 

#band = £ (£k - H )<4(T C k<7 , ( 1 ) 

kcr 

n n 

Hi m p = (So - n) Y, nisa + U £ "«t«;4 , (2) 

i'=l a 1=1 

n 

#hyb = Y,!L Vk ( C ka s to +h - c -) ■ ^ 

i=lka 

c ka creates an electron in the valence p-band at energy £k, n, VCT = s ic s ia is the num- 
ber operator for a Tl s-electron at site i with spin a and energy Eq, U is the (negative) 
on-site interaction, and Vt is the matrix element for the s-p interaction. Its strength 
is characterized by the hybridization function A ( (0) = % Jjj | Vk 1 5 (ffl — Ek a ) ■ We ne- 
glect the k-dependence of Vk setting Vk = Vo but keep the full energy dependence 
of A (a) by using the p-band density of states j¥(ti)) = Y^k^i® ~ £ k) calculated 
from DFT and including relativistic effects J9]- The chemical potential fi determines 
n e = jj Lkcr ( c kcr c kcT) and n s = \ Y%=\ £<? ( n isa)< the average number of p and s elec- 
trons per site. We denote by x = n/N the concentration of Tl impurities. Since Tl 
acts as an acceptor, the ground state corresponds to the Tl 1+ (n s = 2) configuration 
and the Tl 3+ (n s = 0) configuration is split-off from the ground state by the en- 
ergy 8 =£(T1 ) — E(T1 ) > 0. A concentration x of Tl impurities accommodates 
x(n s — 1) electrons (per Tl site), where the number of accepted electrons in the 6s 
level of Tl is measured relative to the neutral Tl 2+ (s l ) configuration having n s = 1. 
These electrons are removed from the valence band leaving behind iiq = 1 — n e 
holes. Thus, the charge neutrality condition reads 



n =x(n s -l), (4) 

which for a given x and temperature T has to be satisfied by adjusting the hole 
chemical potential jj,. Here, we neglect inter-impurity interactions and solve H for 
a finite number of independent negative-^/ centres by using the NRG [37|. For each 
x and each T we need to satisfy (@]i by self-consistently determining the chemical 
potential. In practice we found it more efficient to solve the negative-t/ Anderson 
model on a dense grid of 256 chemical potentials about ji = /i* and to subsequently 
convert these by interpolation to a fixed particle number (i.e. obeying Eq. © for 
each T at given x). 
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3 Transport coefficients 

The transport coefficients for scattering from a dilute concentration x of Tl impu- 
rities are obtained from the Kubo formula 0381 , The electrical resistivity and the 
thermopower are defined by the usual expressions, 

Pimp(r) = -±- , (5) 

e z L n 

S(T) = ~^, (6) 

\e\T L n 

where Ln and L\% are given by the static limits of the current-current and current- 
heat current correlation functions, respectively. In the absence of non-resonant scat- 
tering the vertex corrections vanish and the transport integrals can be written as 

d f^)\ ^„W„ -rW+J-2 



Lij = a deal — J -j—L \^V(cq)t;((q,T)(q' +j -\ (7) 

where <7o is a material-specific constant, f(co) = 1/[1 + exp(co / ksT)] is the Fermi 
function, x{(0,T) is the conduction-electron transport time 11391 

-^—^=2nc imp viA(co,T), (8) 

and A((0,T) =^jImG(fi)i iO + ) is the spectral function of s-holes. The number 
of Tl impurities c; mp per cm 3 is related to x in % by Ci mp = 1 .48 x 10 20 x (using the 
lattice constant «o = 6.46 x 10~ 8 cm of the PbTe rocksalt structure). Eqs.(|7) and d§} 
show that Pi mp (T) depends strongly on the value of A(co, T) around co ~ ji and that 
the sign and the magnitude of S(T) follow from the shape of jy(co)A((0, T) within 
the Fermi window |©| < IkgT . Since PbTe has an unusual non-parabolic band struc- 
ture close to the top the valence band |J9|, and since this region is accessible with 
Tl dopings x < 1%, one sees that the thermopower in particular will be sensitive to 
details of both jY(gj) and A((0,T). 



4 Choice of Model Parameters 

Low temperature Hall effect and tunneling experiments for samples with x > x c 
provide an estimate for jj.* rj 220 me V that we use. Other parameters such as 
Aq = A(n*), required to fix the hybridization function A(co), and U are unknown, 
or they depend strongly on the interpretation of experiments (see discussion be- 
low). However, the measured Kondo-like resistivity for x > x c at low temperatures 
requires that \U\ ~^> -An. We take U/Aq = —8. As argued in 1331 . it is possible 
to estimate U and Aq by interpreting the point contact measurements ll20l within 
our model and by assuming that the Kondo scale T& in the resistivity anomaly 
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at large dopings is comparable to the T c ~ 1.5K, resulting in C/ = — 30meV and 
A = 2.7 meV. While, the precise values of these parameters are required to con- 
sistently explain all measurements, in this paper we shall choose different values 
to those in l33ll in order to substantiate the claim made there, that the qualita- 
tive aspects of those results remain the same for lOmeV < \U\ < 220 me V (and 
U/Aq = —8). Hence, we shall choose Aq = 13.75meV and U = — HOmeV re- 
sulting in 7k ~ 14K, where the Kondo scale T^ is defined via the resistivity 
Pimp(7k) = Pimp{T = 0)/2 for x ^> x c . This is the same order of magnitude as the 
perturbative scale T K = A Q (U /2A ) l / 2 exp(-7tU /8A Q ) = 0.08644 for the negative- 
ly Anderson model I 



5 Qualitative considerations 

Before presenting our results, it is instructive to consider the atomic limit Vq = 0. 
For x = the chemical potential lies in the gap between the valence and conduc- 
tion bands. For finite but very small x each Tl impurity accepts one electron, i.e. 
n s w 2 and hq « x grows linearly with x. At the same time, the chemical poten- 
tial shifts downwards into the valence band jj. < E v , where E v denotes the top of 
the valence band. This implies that the splitting 8(fi) = — (2(et> — fl) + U) between 
donor and acceptor configurations decreases. Eventually, at a critical concentration 
x = x c , the chemical potential reaches jii = n* = Eq + U/2 where 8(h) — 2(/i —y,*) 
vanishes and the system is in a mixed valence state where the Tl 1+ and Tl 3+ con- 
figurations are degenerate. In this situation n s = 1, and any further doping with Tl 
cannot increase the hole carrier density beyond the value no(jU*), i.e. one has self- 
compensation with a pinning of the chemical potential to jj. * [ 32 1 . Experiments show 
that the properties of Pbi_ v Tl v Te change dramatically forx > x c ra 0.3 % but a static 
mixed valence state implied by the above atomic limit cannot capture many of these, 
for example, the anomalous upturn of the resistivity at T < 10K. 

For finite Vo, quantum fluctuations between the degenerate states Tl 1+ and Tl 3+ at 
jj. = ji* become important and lead to a charge Kondo effect ll28l . This significantly 
affects all static and dynamic properties 11291 and needs to be taken into account in 
describing the experiments. This charge Kondo effect is important also for fi > jj.*, 
since a finite charge splitting 8(/i) > in the negative-C/ Anderson model is similar 
to a Zeeman splitting in the conventional spin Kondo effect [41 1. The latter is known 
to drastically influence all properties. Thus, for the whole range of concentrations 
x, one expects dynamic fluctuations to play an important role in the properties of 
Pbi-^TETe. 
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Fig. 1 Hole carrier density no versus Tl doping x in atomic % (at. %) for T = OK and T = 77 K. 
The dot-dashed line is the hypothetical no for one hole per Tl. Filled circles are experimental data 
from Hall number measurements at T = 77 K 1201 . Inset (a): hole chemical potential fl versus x at 
T = OK and T = 77K (dotted line indicates value of the pinned chemical potential fl* = 220meV. 
The value of the charge splitting 8(fl) = 2(fl — fl*) can also be read off from this plot. 



6 Carrier Concentration 



Figure Q] shows the hole carrier density no(x) versus Tl concentration x (in %) at 
T = and at T =77 K. For x < x c ~ 0.3, «o is linear in x, i.e. each Tl contributes 
one hole, as in the case of vacancies on Pb sites (dot-dashed curve in Fig. Q}. For 
x> x c k, 0.3, no(x) saturates rapidly with increasing x for T = and more slowly 
at finite temperature. The hole chemical potential fl, shown in Fig.QJ, grows non- 
linearly with iati< x c and rapidly approaches the value fl* = 200 me V for x > 
x c , both at T = and at T = 77 K. Both of these results are in good qualitative 
agreement with experimental data ll20l (see J33l for comparisons to more data lfT9l ). 
Since we use the realistic band structure, we are also able to obtain quantitative 
agreement for the saturation density «o ~ 0.7 x 10 20 /cm 3 . The self-compensation 
effect at x ^> x c is a "smoking gun" signature for the charge Kondo state: on entering 
this state, the Tl ions fluctuate between Tl 1+ and Tl 3+ so the average valence Tl 2+ is 
neither a donor nor an acceptor and the carrier density cease s to increase for x > x c . 
Counterdoping with donor ions, removes holes, increases fl and the charge splitting 
and consequently destroys the charge Kondo state. Indeed, counterdoping with In, a 
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donor ion, has been shown to destroy the charge Kondo anomalies in the resistivity 
and the onset of superconductivity | 



7 Tl s-Electron Spectral Function 




E-E F [meV] 

Fig. 2 Main panel: 7 = spectral function A(E, T = 0) versus energy E — Ef for a range of Tl 
dopings x. Here, £p = Ep(x) is the self-consistently calculated Fermi level for each x. Inset (a): 
region near E = Ep showing the charge Kondo resonance. Inset (b): valence band local density 
of states ■JV(E') versus E — E-g for x as in the main panel. The conduction band (not shown) lies 
190 meV above the top of the valence band and is not important for the present study. The hole 
carrier densities no for each x, obtained from Fig.fT] are given in the legend. 



Figure |2] shows the evolution with Tl doping of the single-particle spectral func- 
tion A (E, T = 0) at zero temperature. For small doping, the hole chemical potential 
H lies above jJ.*, within the shallow part of the valence band density of states (see 
inset Fig. ^p), and consequently the splitting <5(/i) = 2(/i — ju*) is large. A large 
charge splitting in the negative-i/ Anderson model acts like a large Zeeman split- 
ting in the corresponding positive-t/ model. Consequently for x <C x c , the spectral 
function is strongly polarized, with most weight lying in a Hubbard satellite peak 
far below the Fermi level Ep with no other peaks in the spectral function. With in- 
creasing x, the charge splitting 5(jtt) decreases, the spectral function becomes less 
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polarized, and for x > x c , in addition to the Hubbard peak discussed above, a charge 
Kondo resonance develops close to Ej;. The latter is a result of dynamic valence fluc- 
tuations between the almost degenerate Tl 1+ and Tl 3+ configurations. Point contact 
measurements for PbTe samples doped with x > x c show the existence of two quasi- 
localized states ll20l . a narrow one of width 6meV close to the Fermi level and a 
broader one of width 12meV further below the Fermi level, in broad agreement 
with our theoretical calculations ll33l . Indeed, in [|33l , this interpretation of the ex- 
periments was used to estimate U from the separation of the lower Hubbard band 
from the Kondo resonance. 



8 Resistivity Anomaly 
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Fig. 3 Normalized impurity contribution to the resistivity Pim$(T)/p u versus temperature T and 
a range of Tl concentrations x (in %, see legend). Note that for T > 10K, the total resistivity will 
be dominated by the phonon contribution, not shown in this figure. p„ = 2q mp / (e 1 nKNpOo) is the 
residual resistivity for unitary scatterers with Nf = N(jl") the density of valence band states at 
/! = [l* and o"o a material dependent constant. The vertical arrow indicates the Kondo temperature 
7k = 14K at the charge Kondo degeneracy point. Inset (a): residual resistivity po = xpi mp (T = 
0)/p„ versus x (filled circles). Dashed-line: residual resistivity for static impurities in place of Tl. 



The temperature and doping dependence of the impurity resistivity, pj mp , is 
shown in Fig. [3] For x < x c , Tl impurities act as acceptors with a well-defined va- 
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lence state (Tl 1+ ). They therefore act as weak potential scatterers and consequently 
the resistivity is much below the unitary value, as seen in Fig. [3] Forx > x c , dynamic 
fluctuations between the nearly degenerate Tl 1+ and Tl 3+ states leads to the charge 
Kondo effect and pj mp approaches the resistivity for unitary scatterers at T = 0. For 
x > x c , Pimp is well described by the spin Kondo resistivity 1331 with a logarithmic 
form around T w 7^ ff , where 7£ ff is an effective Kondo scale, and T 2 Fermi liq- 
uid corrections at low T <C 7£ ff , in qualitative agreement with experiment (23J. The 
effective Kondo scale 7£ ff is a function of the charge splitting d(ji) and 7k, and 
approaches the true Kondo scale 7k only asymptotically for x ^> x c (see legend to 
Fig. [3]). Finally, Fig.|3^ shows that the impurity residual resistivity is significant only 
when the charge Kondo effect is operative, i.e. for x > x c , in qualitative agreement 
with experiment [1231 l22l . 



9 Thermopower 

In Figure|4] we show the temperature dependence of the thermopower S(T) for sev- 
eral Tl concentrations in the highly doped regime. At low temperature T <C 7k, the 
thermopower can be obtained from a Sommerfeld expansion of Eq. ©, 

s{T) = -^ t ,r(^-^M), (9) 

3|e| \JY[Ep) A(E F ) J 

where the first term is the band contribution to the thermopower and the second 
term involving the Tl s-electron spectral function is that due to the charge Kondo 
effect of the Tl impurities. Since Tl ions act as acceptors, the Tl ^-electron spectral 
function, A(E), has most of its weight below Ep and its slope at Ep, like that of 
JV(E), is negative (see Fig. [2] and Fig. EJ)). Consequently, at low temperature, the 
charge Kondo and band contributions to the thermopower are both p-type and com- 
pete with each other. This implies that thermopower can undergo a sign change at 
low temperature, depending on details, such as the Tl doping level. Recent measure- 
ments show the occurrence of such sign changes at low Tl dopings [22] consistent 
with a charge Kondo effect. However, our effective single band Anderson model for 
Tl ions in PbTe is not accurate enough to describe the details of these sign changes, 
for the following reason. The valence p-band of PbTe actually consists of two sub- 
bands, light and heavy hole bands due to different hole pockets in the Brillouin zone 
ll42ll . These bands with densities of states JV\{E) and jVi(E) become successively 
occupied with increasing Tl doping, and for x « x c there is significant population 
of both. For simplicity, in Sec. [2] the effect of both bands was approximated only 
in the density of states JV{E) = JY\(E) + JVi(E). For transport properties, such as 
thermopower, which depend sensitively on the relative populations and mobilities of 
such sub-bands, an effective one band model is, in general, not adequate. We expect, 
however, that it is approximately correct in the high doping limit when transport is 
mainly due carriers in the heavy hole band. In this limit, and at temperatures T 3> 7k, 
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when the charge Kondo effect is suppressed, our results for S should reproduce those 
from just the band contribution (,yV(E)) within a constant relaxation time approxi- 
mation x(E,T) = To |0. Fig.|4]shows that this is indeed the case forx > 0.9% and 
T > 175K. The large thermopower in the range 300 — 800K increases with decreas- 
ing carrier density, reflecting the enhancement of jV'(E-p)j jV{E-p) with increasing 
Ep (i.e. decreasing hole chemical potential fi). In Fig. |4j? we show the electronic 
figure of merit zqT = PT/K e , where K e is the electronic thermal conductivity. The 
quantity zqT is an upper bound to the true figure of merit ZT . We find zqT = 2 and 
ZqT = 4 at T = 800K for x ps 1% and x « 2%, respectively, approximately twice 
larger than the reported measured values of ZT at the same temperature and doping 
levels lfT7l. This is reasonable, given our neglect of the lattice thermal conductivity. 
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Fig. 4 Main panel: Thermopower S versus temperature at different Tl concentrations x (in %, see 
legend) in the highly doped regime x > x c . The two solid lines labeled by «o correspond to S for 
the two highest doped cases (i.e. for x = 0.9 and x= 1.8) assuming a constant transport time. 
Inset (a) shows the temperature dependence of the dimensionless electronic figure of merit zoT for 
x > x c = 0.3 (in %, see legends). This is indicative of trends only, since the lattice contribution 
to the thermal conductivity has been neglected in zoT. The vertical arrow denotes the particular 
choice of Kondo temperature 7k = 14 K made in this paper (see discussion in Sec.[4}- 



Within our model calculations, the charge Kondo contribution to the thermopower 
at low temperatures competes with the main contribution coming from the band 
instead of supplementing it. However, the contribution from the charge Kondo ef- 
fect is generically quite large, and can reach values of order 76/X V/K for splittings 
5(ju)/7k > 1. This is indicated in Fig.|5]for the charge Kondo thermopower of the 
negative-^/ Anderson impurity model ||29l using a flat band appropriate for impuri- 
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ties in metallic systems. Hence, in principle, impurities exhibiting the charge Kondo 
effect could enhance the overall thermopower by the above value. This could offer 
a route to further improving the thermoelectric properties of appropriate materials. 
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Fig. 5 Thermopower S versus reduced temperature T/Tg_ of the negative-?/ Anderson impurity 
model in a metallic host |29l for \U\/Aq = 8 and several values of the charge splitting 8(fl). The 
latter is denned as the splitting between the empty and doubly occupied states of the negative- 
U Anderson impurity model and is analogous to the self-consistently calculated S(fl) for the Tl 
impurities in PbTe discussed in this paper. 7k denotes the Kondo scale. The charge Kondo ef- 
fect provides a mechanism for large thermopowers of order k^fe w 76/1 V/K at charge splittings 
5(jU)/7k > 1. Since Tk can be made small, due to its exponential dependence on U and Aq, such 
splittings are accessible in potential realizations of the charge Kondo effect, e.g. in negative-?/ 
centers in semiconductors or in molecular junctions 11291 . 



10 Conclusions 



In summary, we investigated the low temperature properties of Pbi_. r Tl r Te within a 
model of Tl impurities acting as negative U centers. Our NRG calculations explain 
a number of low temperature anomalies of Pbi_ r TLTe, including the qualitatively 
different behavior below and above the critical concentration x c , where x c ~ 0.3 %. 
They support the suggestion that the charge Kondo effect is realized in Pbi_ v TLTe 
11231 [32ll . At x = x c , two nonmagnetic valence states of Tl become almost degen- 
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erate and the ensuing pseudospin charge Kondo effect results in a Kondo anomaly 
in the resistivity for x > x c and a residual resistivity approximately linear in x. Our 
results for these quantities and the carrier density hq (x) are in good qualitative agree- 
ment with experiments J23H T9 22 201. For the Tl ^-electron spectral function, we 
predict that one peak should be present far below Ep for x < x c and that a second 
temperature dependent Kondo resonance peak develops close to, but below Ep, on 
increasing x above x c . This provides a new interpretation of measured tunneling 
spectra ll20l . which could be tested by temperature dependent studies of tunneling or 
photoemission spectra. We also showed that the competing charge Kondo and band 
contributions to the low temperature thermopower imply that there are sign changes 
in the thermopower at low Tl concentrations. Investigating these in detail, particu- 
larly at x < 0.3%, could shed further light on the charge Kondo effect in this system. 
In the future, it would be interesting to extend this work to investigate the effects of 
disorder, phonons, light and heavy hole bands, and non-resonant scattering channels 
on the thermoelectric properties of Tl doped PbTe. 
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